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Introduction 

Solar cells made from gallium arsenide (GaAs), with a room temperature bandgap 
of Eg = 1.43 eV have exhibited the best measured open circuit voltage (V oc ) of 1-05 
V at 1 AMO, 25C [ref. 1]. The material InP is in many ways similar to GaAs. A 
simple calculation comparing InP to GaAs then shows that solar cells made from 
InP, with Eg = 1.35 at 300K, should exhibit the best measured V oc of ~950 mV at 
1 AMO, 300K. However, to date, the best measured V oc for InP solar cells made by 
any fabrication method is 899 mV at AM1.5, 25°C [ref. 2], which would translate to 
912 mV at 1 AMO, 25°C. 

The Voc of an n+pp+ InP solar cell is governed by several factors. Of these, some 
factors, such as the thickness and doping of the emitter and base regions, are easdy 
controlled and can be adjusted to desired values dictated by a good performance- 
optimizing model. Such factors have not been considered in this investigation. There 
are other factors which also govern V 0 c, and their values are not so easily controlled. 
The primary ones among these are 1) the indirect or Hall-Shockley-Read lifetimes in 
the various regions of the cell, 2) the low-doping intrinsic carrier concentration n^ of 
the InP material, 3) the heavy doping factors in the emitter and BSF regions, and 
4) the front surface recombination velocity Sp. We have investigated the influence of 
these latter factors on the V 0 c of the n + pp+ InP solar cell and have used the results 
to produce a near-optimum design of the n + pp + InP solar cell. 

Theoretical Approach 

We have developed a fairly comprehensive, one dimensional, closed form solution, 
computerized model of the shallow homojunction InP space solar cell. This model 
assumes a cell with three homogeneously doped regions, an n + emitter, a p-type base 
and a p+ BSF region. Details of the model have been presented elsewhere [ref. 3]. In 
our earlier near-optimum design of the n + pp + shallow-homojunction InP solar cell 
using this model [refs. 3,4] no accounting was made of the heavy-doping effects in the 
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heavily doped emitter and BSF regions. As of this day, to the best of our knowledge, 
measured heavy doping factors or equivalent bandgap narrowing values as functions 
of doping for n-type and p-type InP have not been published in the literature, as 
they have been for GaAs [ref. 5]. This lack of information makes it difficult to take 
into account heavy doping effects in InP cell modeling. Even so, it is now realized 
that it is very important to take these into account in order for the modeling to be of 
much value. Hence, some way must be found to reasonably estimate the heavy doping 
effects. A prior attempt [ref. 6] at this has been to represent the heavy doping effects, 
including Fermi- Dirac statistics, by an equivalent effective bandgap narrowing AE g 
in each heavily doped cell region (emitter, BSF) and to estimate the amount of AE g 
in these regions by matching the calculated and measured internal spectral response 
in the short and long wavelength region of the incident light spectrum. We feel that a 
less roundabout and more accurate method is to lump all the heavy doping effects into 
an effective intrinsic carrier concentration nj e and to thereby define a heavy doping 
factor HD in each heavily doped region as, for example in the emitter, by 

HD e = ^ 

n i0 

where n ieE is the the effective value of n, in the emitter and n i0 is the low-doping 
value of n; in InP. A similar heavy doping factor HD B g E 1S defined for the BSF region. 
The values of HD E and HD B sf are then estimated by matching the calculated and 
measured curves of not only the internal spectral response but also the illuminated 
I-V and log Isc “ Voo 

The above approach still leaves one major problem. It allows the estimation of 
IID e and HD B g E , eac h for only one value of doping, that of the experimental cell 
whose measured curves w r ere matched. However, what is needed for general modeling 
purposes is HD E and HD BS f as functions of their respective dopings: in other words, 
curves of HD| and HD| gF vs. their respective dopings, as are available for GaAs at 
300K [ref. 5], as shown in figure 1. In the absence of available measured data on InP, 
we made the reasonable assumption, as a first approximation, that for both n-type 
and p-type InP, the nature of the HD 2 versus doping curve is similar to that of GaAs. 
With this assumption, the curve for InP can be obtained from that of GaAs if one 
point on the curve could be determined. As stated earlier, this can indeed be done 
by matching the calculated and measured curves of illuminated TV, log I sc - V 0 c and 
internal spectral response. 

Our computer model was designed from the start to correctly take into account 
heavy doping effects in all three regions of the cell, if necessary, by means of appro- 
priate boundary conditions involving the heavy doping factors in each region. The 
boundary condition at the emitter-base junction is that of the generalized ‘law of the 
junction’, valid at all injection levels, while that at the base-BSF interface is that of 
the low-injection effective interface recombination velocity S pp+ given by 
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[1] 


s = Na > base / HD bsf \2 P n ,BSF 

pp+ N a B SF VHD base / L BSF tanh( L WfiSF ) 

where N a is the doping, HD is the heavy doping factor, as defined earlier (e.g. HD B g F 
= n ie,BSF/ n io)> an( i D, L, w are the minority carrier diffusivity, diffusion length and 
width of the region, respectively. For a base doping of less than 10 17 cm -3 , HD base = 
1, and it is seen in equation [1] that the net effect of heavy doping in the BSF region 
is to increase the effective base- BSF interface recombination velocity S pp + by a factor 
of (HD bsf ) 2 . As we shall see from the calculated results in the following section, the 
overall effect of heavy doping in the BSF region is to increase the minority carrier 
recombination in the base at the base-BSF interface, thereby increasing the dark 
saturation current in the base and reducing the open circuit voltage of the cell. 

As to the task at hand, that of finding HDg and HD B g F at one doping, we 
chose to match theoretical and measured curves of illuminated I-V, log I sc - V oc and 
internal spectral response for a 17.9% efficient n+pp+ solar cell labeled Spire 6 [ref. 
7], made by the Spire Corporation, for which such measured curves were available. 
All measurements on this cell were made at the NASA Lewis Research Center in 
Cleveland, Ohio, and were made available to us by courtesy of Mr. Russell E. Hart. 

We obtained excellent matches between the calculated and measured curves of 
Spire 6 for specific values of IID E = 1.4 for the n-type emitter doping of 1 x 10 18 cm -3 
and HD bsf = 4.87 for the p-type BSF doping of 5 x 10' 8 cm -3 . By noting that each 
of the two HD 2 versus doping curves had to pass through its corresponding specific 
point indicated above and also that each HD has to become unity at dopings of 10' 7 
cm -3 and below, we were able to express the general variations of IIDg and HD B g F 
versus their respective dopings in equation form as: 


HD E,n-lnP — !-92 [HD 2 _ GaAs - 1] + 1 [2] 

except when HD n _ GaAs drops below unity, in which case, for simplicity, we used the 
approximation 


HD|,„. InP = hd;_ G iA , 


[3] 


Also, 


HDbsf, p-inP - 6.1 • [HD 2 _ GaAs - 1] + 1 [4] 

The +1 on the right-hand side of equations [2] and [4] ensures that HD in InP becomes 
unity when IID in GaAs becomes unity, at the lower dopings. 
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The above equations were incorporated into the computer model so that, for any 
doping concentration in the emitter and BSF, the corresponding heavy doping factor 
could easily be calculated. This capability to calculate the heavy doping factors 
HD e , HD bsf as functions of doping was the only new addition that our computer 
model required. The computer model was already capable of calculating all other 
geometrical and material parameters needed to simulate the solar cell, as in the past 
[refs. 3,4]. 


Calculated Results and Discussion 

In order to see how the V oc varied with each of the individual parameters expected 
to affect it, we first chose a somewhat modified version of the Spire 6 cell design to 
serve as a baseline cell. The primary modifications were: a) a reduced grid shadowing 
from 4.8% of Spire 6 to 4.0% of the baseline design, b) a thinner emitter, from 400A of 
Spire 6 to 200A of the baseline cell, c) a heavier emitter doping, 5 x 10 18 cm -3 from 
1 x 10 18 cm -3 of Spire 6, d) a thinner base, 2 /im instead of 3 /xm, and e) a lighter BSF 
doping of 2 x 10 18 cm" 3 for the baseline cell instead of 5 x 10 18 cm" 3 for Spire 6. These 
modifications reflected current practice in cell design and resulted in an improved cell 
efficiency of 19.23% from the 17.9% of Spire 6. The baseline design also resulted in 
a somewhat larger V 0 c of 873 mV from the 868 mV of Spire 6. Table 1 shows all the 
geometrical and material parameters of the Spire 6 and the baseline solar cells. Note 
that the numbers in the baseline cell column are for 300 K as opposed to 298 I\ for 
Spire 6. This results in a significantly longer radiative lifetime in the base for the 
baseline cell, since the base radiative lifetime is extremely sensitive to temperature. 

In our model, the indirect or Hall Shockley-Read (HSR) lifetime in any cell region 
is assumed to be inversely proportional to the doping in that region. This assumption 
is tantamount to the assumption that indirect recombination in InP occurs via a single 
dominant recombination level whose density is proportional to doping density. Thus, 
in any cell region, 


t HSR - 


Fhsr 

N 


[5] 


where r HS R is the constant of proportionality between t HS r and the reciprocal of the 
doping concentration N. Only two different values of TjjgR are needed, one for n-type 
and one for p-type InP. These were obtained from our matching the calculated and 
measured curves for Spire 6, as were the values of HD E , HD BSF and^S F . For the 
baseline value of n ; , we chose the prevalent value of ^(300 I<) = 1.2 x 10'cm -3 , which 
may be on the high side according to current thinking [ref. 2]. 

Thus, knowing the baseline values of HD E , HD B gFi Thsr,p) Fhsr.ih S f and n ; , 
all except n; being derived from the Spire 6 match, we could calculate the complete 
performance of the baseline cell at 300 K. The second column of Table 1 gives the 
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geometrical, material and performance parameters of the baseline cell. Note that the 
baseline cell has a thinner, more heavily doped emitter, a thinner base and a less 
heavily doped BSF layer than the Spire 6 cell. 

Next, we made a large number of parametric variation runs of our computer 
model simulating the baseline cell, varying one parameter at a time from among 
IIDg, HD bsf , r B g B p and r B g B n (in the form of relative HSR or indirect lifetime 
Tj/Ti 0 in the n-type emitter and p-type base and BSF regions), the low-doping intrinsic 
carrier concentration nj 0 and finally the front surface recombination velocity Sp. It 
should be noted that in our model, Sp is the area-weighted average effective surface 
recombination velocity over both the part of the front surface which is not in actual 
physical contact with the front grid metallization and the part which is in such contact 
[ref. 3]. 

The results of the parameter variation runs are shown in Figs. 2,3,4 and 5, which 
show, respectively the variation of V oc versus 1) the HSR relative lifetime Tj/r i0 in 
the emitter, base and BSF, 2) the low-doping intrinsic carrier concentration n^, 3) 
the BSF heavy doping factor IIDggp and 4) the effective front surface recombination 
velocity Sp. 

We now discuss each of these figures in turn. In Fig. 2, we see that if, by using 
better material and device processing, we can make the indirect or HSR lifetimes in 
all regions of the cell to be ten times longer than their current-technology values in 
Spire 6, then V oc can be improved by slightly over 10 mV, which, while not being a 
very large gain in V oc still is a significant gain. Fig. 2 also shows that V oc begins 
to level off as the HSR lifetime is made even longer. This is because when the HSR 
lifetime becomes very long, the overall lifetime is controlled by the radiative lifetime 
which, for a given doping and temperature, is an inherent property of the material InP 
and is independent of the quality of the starting material or of the device processing 
parameters. Thus, it appears from Fig. 2 that by making the HSR lifetime in each 
region indefinitely long by using the highest quality starting material and using the 
cleanest possible device processing, the most to be gained in V oc is about 15 mV or 
so. 


Next considering Fig. 3, we see the variation in V oc as the 300 K, low-doping 
intrinsic carrier concentration n; 0 decreases from 1.2x 10 7 cm -3 to a hypothetical value 
of 5 x 10 6 cm -3 , to show the nature and extent of variation of V oc with nj 0 . If n;o 
were somehow halved from its presently used value of 1.2 x 10 7 cm -3 to 6 x 10 6 cm -3 , 
then Fig. 3 shows that V oc would increase by ~25 mV. The values of n; 0 and V oc 
indicated by the two arrows in Fig. 3 are as follows. The one corresponding to njo 
= 1.2 x 10 7 cm 3 represents the prevalent value of n; based on both the density of 
states functions N c , N v , and the bandgap E g corresponding to InP. The other point 
is the value of nj 0 that would be obtained by using the density of states values of 
GaAs but the bandgap of InP-in other words, the value of n; that InP would have if 
its N C N V product were the same as that for GaAs. The difference in V oc between the 
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two points is about 18 mV, and it indicates that, in comparison with GaAs, InP has 
an 18 mV V 0 c disadvantage merely because its N C N V product is too large compared 
to that of GaAs. 

There is some uncertainty about the value of nj in InP. Yahia and Coutts have 
done an extensive investigation of this problem, and believe that n; in InP at 300 
K should be around 8 to 9 x 10 6 cm' 3 [ref. 8]. If this is true, then the expected 
maximum V 0 c of InP would be higher by about 10 mV. 

In Fig. 4, we show the calculated variation of V 0 c for the baseline InP cell with 
the BSF heavy doping factor HD E sf = n ie,BSF/ n iO The point corresponding to the 
heavy doping factor calculated using equation (4) is indicated by the arrow. It is seen 
that the lower the heavy doping factor, the higher the V oc - 

If the qualitative behavior of heavy doping factor versus doping in p-type and n- 
type InP is indeed similar to that for GaAs, as we have assumed, then it would appear 
at first sight that it would be advantageous to have an n-based cell with n + BSF. 
However, only detailed calculations would show whether that would be the case or 
not, for the following reasons: First, in a p+nn"*" InP cell, the p + emitter would have 
to be relatively thick in order for the low hole mobility in it to not cause an excessively 
large emitter sheet resistance component of series resistance. With a thicker emitter, 
a significant portion of the total cell photocurrent would now come from the emitter, 
as would the dark saturation current. In that case, any gain due to a low (lower than 
one) HDbsf could be more than counteracted by the high HD E in the p + emitter. 
We are in the process of doing such detailed calculations for the p^nn" 3 " cell. 

Finally, Fig. 5 shows the variation of V oc versus the effective front surface 
recombination velocity Sp. Here, it is seen that over its entire range of values 
from Sp < D pE /L pE , the diffusion velocity of minority carriers in the emitter, to 
Sp >• D pE /L pE , Voc degrades by only about 3 mV. The reason why V 0 c is relatively 
insensitive to Sp is that in any relatively well-designed InP solar cell, as in the base- 
line cell, the V oc is controlled primarily by the dark saturation current J 0 in the base, 
which is several times (~10 or more) larger than the Jo from the emitter. Thus, 
the Sp variation, which may affect Jo, emitter by a f ac t° r °f a f ew i does not affect the 
overall J 0 , and therefore V oc , significantly. It is to be noted, however, that Sp has a 
significant effect on the short circuit current Jsc and on the efficiency, therefore, it is 
most desirable to reduce Sp to a very low value (Sp < D pE /L pE ) by means of front 
surface passivation or the use of an appropriate higher bandgap window layer. 


Near— Optimum Cell Design 

Having investigated the individual variation of V oc with each of the parameters 
F HSR.n an d Thsr.pi n i0 , IID bsf and S F and some other variations (V oc vs. base 
thickness, base doping, emitter thickness, emitter doping and HD E ) not shown here, 
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we next proceeded to design a near-optimum InP n+pp+ solar cell for 1 AMO, 300 K 
operation with a view to maximizing its efficiency (not V oc ) 

Our approach to the near-optimum design was to ignore limitations imposed 
by current technology on such parameters as r HSR p , r HSRn and Sp. With improved 
technology, in the near future, r HSR p and r HSR n could realistically be about ten times 
the values obtained in the Spire 6 match, while Sp could very likely be ~ 10 3 4 cm/s or 
less. Thus we have used values for Thsr.p an d Thsr,,, that are ten times the values 
obtained in the Spire 6 match and used in the baseline cell. This is tantamount to 
assuming that in the near future, material and processing quality is likely to improve 
to the extent that at moderate and high dopings, the overall lifetime in InP will be 
limited by the radiative lifetime. We have also used Sp = 10 4 * cm/s, again assuming 
that the significant amount of research on surface passivation and searching for an 
appropriate window material will result, in the near future, in an effective front SRV 
of S F = 10 4 cm/s or less. As to the technology- independent parameters, such as 
HDp, HD R sp and n; 0 we have used the same values as used earlier, namely, HDp and 
HDggp based on equations [2], [3] and [4], and n i0 = 1.2 x 10 7 cm -3 at 300 K. 

Table 2 gives the complete geometrical, material and performance parameters of 
the near optimum cell. Note that this optimized cell has a realistic expected efficiency 
of 22.61% and a V oc of 901.6 mV. According to our model, higher V oc values are 
obtainable at some sacrifice in J sc and efficiency. 

Because of heavy doping effects, our near-optimum design has relatively lighter 
dopings in the BSF and base. It may seem at first sight that with a BSF doping only 
ten times the base doping, the base/BSF interface would not be an effective reflector 
of minority carriers. However, higher BSF dopings result in the base/BSF interface 
being an even worse reflector of minority carriers because of heavy doping effects in 
the BSF. 


Concluding Remarks 

1. The primary factor limiting the V oc in an n+pp+ solar cell is the heavy-doping 
factor in the p+ BSF. 

2. In a well-designed n + pp + InP solar cell, the effective front SRV Sp has a 
significant effect on I sc and T) but a relatively minor effect on V oc - 

3. If the effective front SRV is reduced to ~ 10 4 cm/s and the HSR lifetime 

coefficients r R g R p , r R g R n are increased by a factor of 10, then it is theoretically 

possible to obtain a V oc of 902 mV and T) = 22.6% at 1 AMO, 300 I\, for an n+pp+ 

InP solar cell. ENen higher V oc values should be obtainable at a sacrifice in I sc and r/. 
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